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Abstract 

O _ 

CN ■ We reconsider chiral perturbation theory in a finite volume and develop a new computational 

^ I scheme which smoothly interpolates the conventional e and p regimes. The counting rule is kept 

' essentially the same as in the p expansion. The zero- momentum modes of Nambu-Goldstone bosons 

are, however, treated separately and partly integrated out to all orders as in the e expansion. In this 
new scheme, the theory remains infra-red finite even in the chiral limit, while the chiral-logarithmic 
effects are kept present. We calculate the two-point function in the pseudoscalar channel and show 
that the correlator has a constant contribution in addition to the conventional cosh function of 

m 

^> . time t. This constant term rapidly disappears in the p regime but it is indispensable for a smooth 

I convergence of the formula to the e regime result. Our calculation is useful to precisely estimate 

the finite volume effects in lattice QCD simulations on the pion mass M.^ and kaon mass Mk, as 

in 

O ■ well as their decay constants F^^ and Fk- 



X 
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I. INTRODUCTION 



Recent progress in lattice QCD has made it possible to simulate QCD in a realistic set- 
up, i.e. with the (2+l)-fiavor sea quark masses near the physical point. As the precision of 
the data analysis goes high, however, more precise study of systematic effects is required. 
Finite volume effects are particularly important when quark masses are reduced to near the 
chiral limit, since the correlation length of the system rapidly grows, which is induced by 
the dynamical chiral symmetry breaking 

The chiral symmetry breaking makes a mass gap between the Nambu-Goldstone bosons, 
which eventually become massless in the chiral limit, and the other hadrons, which retain a 
mass around the QCD scale Aqcd- K is, therefore, the pions that are the most responsible 
for the effects of the finite volume V when the size of the system L or V^^^ is well above 
I/Aqcd- 

With this motivation, a number of studies have been devoted to understand the finite 
volume effects within the theory of pions, which is known as chiral perturbation theory 
(ChPT) {2], Using the lattice data for the low-energy constants as inputs, one can quantify 
the finite volume effects from the pion fields. These studies are also useful for improving the 
determination of the input low-energy constants themselves. 

To investigate ChPT in a finite volume, two perturbative approaches have been proposed 
so far. One is the p expansion j4-7|, which has just the same form as the perturbative series 
in an infinite volume, but momentum integration is performed in a discrete space in the 
units of 1/L. Denoting the mass of a generic (pseudo) Nambu-Goldstone boson by M, this 
p expansion is valid when ML ^ 1, which is called the p regime. 

A nonperturbative technique is required when ML ^ 1 (the e regime) since the zero- 
mode's contribution to the propagator of the pseudo Nambu-Goldstone bosons blows up and 
fluctuation ~ cannot be perturbatively treated, which is well-known as the critical 

fluctuation due to the symmetry breaking. A solution to this problem was given in terms of 
ihe so-called e expansion in Refs. js iSl and later the study is extended in various directions 
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25| . In this scheme the zero-momentum mode is separately treated and integrated out 



exactly, while all the remaining non-zero momentum modes are treated perturbatively. Since 
the e expansion treats the mass term as a next-to-leading order (NLO) contribution, the 
number of terms in the chiral Lagrangian is reduced compared to the p expansion and the 
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typical chiral-logs are invisible in the calculation at NLO. Note here that the exaci integration 
here refers to the term that is leading order in the quark masses m. 

One may ask what happens in between: when ML ~ 1. The answer should be given in 
either ways of the expansions since the p and e expansions should eventually converge to 
give the same result as the order of loop expansion increases. But it is difficult already at 
the two- loop level, to confirm such a convergence between the p regime [3| and e regime 25 1 
calculations unless one directly checks the numerical values, since their analytic forms look 
quite different. It is, therefore, important and useful for the practical calculation, to find 
a new way of expansion which smoothly interpolates the p and e expansions while keeping 
the calculation at the one-loop level. Intuitively, this one-loop level interpolation should be 
possible in the simplest way, by keeping all the terms that appear in the NLO Lagrangian 
in both expansions. 

In fact, such a calculation is demanding. Although recent developments in computational 
facilities have allowed us to simulate unquenched lattice QCD near the chiral limit, it is still 



difficult to fully satisfy the condition ML ^ 1. On t 



le other hand, no study has until now 



reached deep inside the e regime keeping ML <g 1 26 



33| . Although results have often been 



compared favorably to the e expansion of ChPT, there may still be large systematic errors 
due to the condition ML <^ 1 not being well fulfilled. 

Recently a new approach which smoothly connects the p expansion and e expansion (and 
which remains valid even in the region ML ~ 1) was proposed in Ref. 3J]. The new 
prescription is to keep the counting rule of the p expansion but treat the zero-mode non- 
perturbatively as in the e expansion. This new expansion was applied to the calculation of 
the chiral condensate (and the spectral density of the Dirac operator) to NLO and successful 
in maintaining the features of the both regimes: non-perturbative behavior of the zero-modes 
and chiral logarithms. The results are kept infra-red (IR) finite even in the chiral limit 35- 

n nn 

|37| and show a good convergence to the conventional result [38|, l39| in the p expansion for 



the large (valence) quark mass region. A good agreement with a lattice QCD calculation 
was reported in Refs. 40|, l41 |. 



In this paper, we extend the calculation of Ref. [3J] to the two-point functions in the 
pseudoscalar channel. We find that the correlator is expressed by a simple hyperbolic cosine 
function of time t plus an additional constant term, which smoothly connects the conven- 
tional p regime results and those in the e regime. The constant contribution is a peculiar 
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feature of the e expansion. We find that this constant is indispensable to keep the correlator 
IR finite, and show how and where it becomes negligible as entering the p expansion regime. 
Our results are useful to precisely estimate the finite volume effects in lattice QCD on the 
pion mass and kaon mass M^, as well as their decay constants and F^- 

The rest of our paper is organized as follows. In Section HIl we describe in detail our 
new perturbative counting rule in ChPT and the computation scheme which consists of 
three steps. For the first step, the chiral Lagrangian in terms of non-self-contracting (NSC) 
vertices (whose definition is given in the following sections) of non-zero momentum modes 
is calculated in Section lllll The second step is to collect the one-loop diagrams of the 
correlator and perform the non-zero mode's perturbative integrals (Section llVp . The final 
step is non-perturbative zero-mode's integration in Section |Vl The results for the two-point 
functions in the theory with a general number of flavors are presented in Section I VI I (see 
Eq.( pO|) ). For more practical uses, explicit formulas for the Nf = 2 and 2+1 cases are given 
in Section FVIII (see Eq. (lll2p ) as well as how to compare the results with the lattice QCD 
data. Our calculation suggests that there exists a simplified short-cut prescription which 
reproduces the same results. We discuss this simplified scheme in IVIIII Conclusions are 
given in Section HXl 



II. NEW CHIRAL EXPANSION AT FINITE VOLUME 



In this section we review the new counting rule of chiral perturbation which was first 



proposed by Ref. 3J]. We also present our strategy for the calculation of two-point functions. 
We consider an A^/-fiavor chiral Lagrangian in a finite volume {V = L^T), 

C = ^Ti[d^U{xyd^U{x)] - -Tr[Mte^9/7V/[/(x) + UixYe-'^/^f M] + ■ ■ ■ , (1) 

where U{x) G SU{Nf) and 6 denotes the vacuum angle, while S is the chiral condensate 
and F denotes the pion decay constant both in the chiral limit. We note that the higher 
order terms are not explicitly shown here but exist, which is indicated by ellipses. 

In the partially quenched case, we use the replica method where the calculations are 



done within an {Nf + + {N — A'^,))-fiavor theory and the limit — is taken j42|-|44| \ 



^ We do not consider the fully quenched theory in this work. We thus have Nr > in all that follows. 
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Physical unquenched A^/-flavor theory results can be obtained by simply taking = "^f 
where m/ is one of the physical quark masses. 

For the mass matrix, we thus consider a general non-degenerate form: 

M = diag(m^^ , ■ ■ ■ m^^ ■ ■ -, m„,mrf , mg, ■ ■ ■ ). (2) 

Ni N2 Nf 

where we have N = Ni + N2 replica flavors and A'^^ physical flavors. Since our target is a 
single meson system which consists of two quarks, we have written the valence part as if 
there were two different sets of degenerate flavors, where each of Ni quarks have a degenerate 
mass my-. For each valence flavor, the Ni ^ limit has to be taken in the end of calculation 
to complete the partial quenching. 



We parametrize the chiral field in the same way as the e expansion js], by factorizing it 
into the zero- momentum mode Uq and non-zero modes ^{x), 

U{x) = UoexpiiV2^{x)/F). (3) 

In our calculation, we perform exact group integration over Uq, while ^{x) is perturbatively 
treated always imposing 

J d'x ^{x) = 0, (4) 

to avoid double counting of the zero-mode. 

It is known that group integration over U{Nf) manifold is easier and can be analytically 
expressed in a simpler form than the SU{Nf) group case. For this practical reason, we 
consider sectors of fixed topology Q, which is obtained by the Fourier transform of the 
partition function, 

l-j^JdOe'^'^jvUe-''. (5) 

We then absorb the 9 integral to the zero-mode sector: e^^^^fUo — )■ Uq and extend our 
integration to U{Nf) (or U{Nf + A^) in the partially quenched case) group. The phase 
factor in the Fourier transform becomes e*^*^ = (det Uq)^. The conventional = vacuum 
result is obtained by summing each topological sector with a weight given by the partition 
function, which will be discussed later in Section I VII 

We give the same counting rule as in the p expansion for the ^ fields and other parameters, 

d,r^O{p), ax)r^O{p), Mr^Oip'), T , L ^ O {I / p) , (6) 



in units of the cut-off AnF. We assume as usual that the hnear sizes of the 4-dimensional 
volume, L and T, are much larger than the inverse QCD scale Aq^j-, so that the effective 
theory is valid. 

According to the counting rule Eq. let us expand the Lagrangian 



1„ ..o S. 



1 - r2rt2i 
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2 V 2 

t 

+^^Tt[M\Uo - l)e + eiU^o - l)M] + ■ ■ ■ , (7) 

where M^j = {jrii + ■mj)T,/F^. Here we have separated the mass term into three pieces. 
The first one (the second term) gives a non-perturbative weight in the zero-mode path 
integration as in the e expansion and the second one (the third term) has the same form as 
the conventional mass term (of ^) in the p expansion. 

The last term in Eq. ([7]) is a mixing term between the zero and non-zero modes, which 
is unfamiliar either in the e and p expansions. In fact, this term plays a crucial role in 
connecting the e and p regimes. We can treat this mixing term as a perturbation: it is not 
difficult to check 

MiUo-l)r^Oip'), (8) 

and, in particular, a Hermitian combination 

M{Uo + U^o-'^)-0{p'), (9) 

hold in both of the e and p regimes. For some specific cases, by a direct group integration, one 
can confirm that these countings are kept even in the intermediate region where MijL ~ 1 



34l |. We therefore treat Eq. (jH]) and ([9]) as the additional counting rules and treat the last 
term in Eq. ((Tj) as an 0{p^) contribution. These additional counting rules Eq. (jH]) and 
are also supported by the equipartition theorem of energy, where the potential energies of 
weekly interacting system are uniformly and therefore, mass-independently distributed. 

In Table [H we summarize the difference of the three e, p, and our new i (=interpolating) 
expansions of ChPT. 

In the following sections, we calculate two-point correlation function of the peudoscalar 
operators in three steps. For the first step (Section IIIII) . we rewrite the chiral Lagrangian 
in terms of non-self-contracting (NSC) vertices of ^ fields. This corresponds to partly per- 
forming one-loop integrals in the vertices in advance. By doing this, one can renormalize 
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expansion 


parametrization 


counting rule 


e expansion 


Uix) = Uoexp (i^) 


Uo ~ C ~ 0{1/L), M ~ 0{1/L^) 


p expansion 


C/(x)=exp [i^) 


^ ~ 0{1/L), M ~ 0(1/L2) 


i expansion 


C/(x) = ^7oexp (i^) 


[/o ~ C ~ 0{1/L), M ~ 0(1/^2), 
M{Uo - 1) ~ 7W(C/o + Ul - 2) ~ 



TABLE I: Three expansions of ChPT at finite volume. The counting rules are compared in the 
units of the smallest non-zero momentum l/L. Our new expansion in this paper is denoted by "i 
expansion" . 

the coupling constants and the wave function at NLO before starting the comphcated cal- 
culation. Then the second step for the two-point functions (Section IIVI) becomes clearer: to 
collect the remaining diagrams, namely those without self-contractions in vertices, which is 
expressed by the already renormalized quantities, and perform ^ integrals. The third and 
final step is to perform nonperturbative Uq integrals. 

For the perturbative calculation of fields, we use the same Feynman propagator as in 
the p expansion except that the zero-momentum mode contribution is removed: 



(10) 



where (■ ■ ■)g means an integral over ^, whose general expression will be discussed later in 
Sec. IIVI Note that the second term comes from the constraint Tr ,^ = 0. The propagators 
A and G are given by 



1 ^ e-' 



G{x,MiMi) = -y: 



(11) 

(12) 



where the summation is taken over the non-zero 4-momenta 

p = 27r{nt/T,n^/ L,ny/L,njL), 



(13) 



with integer n^'s except for p = (0,0,0,0). For the following calculations, where a non- 
degenerate set of valence and sea quark masses is taken, it is convenient to define a quantity 

1 



A{x,MlM ^ ^ G{x,MlM 



33 > 



G{x,Mt,,Mtd + G{x,M^^,M^^, 



(14) 



Note that both A{x, Mf^, M^ ) and its second derivative 



dlA{x,MlM: 



3J> 



Mm{x,MlM: 



J0> 



(15) 



are UV finite even in the hmit x = 0. Also, note that both vanish when M?- = Mf,. 

As a final remark of this section, we note that the above parametrization Eq. Q gives 
rise to a non-trivial Jacobian in the functional integral measure. It is uniquely determined 
by the left-right invariance of the group integrals. A perturbative calculation 
shown that the Jacobian is expressed by 

Nf 
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17| has 



J{U,,i) = exp(- j d'x^^Tie 



(16) 



to 0{p^). It plays a role just as an additional mass term in our calculation. 



III. CHIRAL LAGRANGIAN AT ONE-LOOP 



Since our target system is a complicated mixture of Uq matrix model and perturbative ^- 
fields, we first simplify the chiral Lagrangian and collect relevant pieces for our computation. 
In particular, by introducing non-self-contracting vertices, we can renormalize (at the one- 
loop level) the coupling constants and the ^ fields in advance. 



A. Next-to-leading order (NLO) terms 



Without source terms, we have eight NLO terms, whose low-energy constants are denoted 
by Lj's (i = 1, ■ ■ - 8) In our perturbative expansion at 0{p^) and 0{p^), the terms with 
Li, L2, L3 (and the Wess-Zunimo-Witten term |45|, |46[ as well) do not contribute to pseudo- 
scalar meson masses and decay constants. By explicitly expanding U{x) = Uoe'^^^^^^^^ in 
^, it is sufficient to consider 



Cnlo = --Tt[M^Uo + U^M] X 



+E 
+E 



^jr,[UY.M'ff + L,Ml 



f 



X 



16 
F2 



LeY^Mfj + LsMl 
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Note that we can always omit the constant terms unless source terms are inserted (the source 
insertion is separately discussed below). It is also important to note in the above expansion 
that the only Lq term has non-trivial Uq dependence at 0{p^). 



B. Non-self-contracting (NSC) vertices 

For one-loop level calculations, it is convenient to rewrite the chiral Lagrangian so that 
quantum corrections are partly included. This is performed by simply adding and subtract- 
ing all possible ^-contractions of the n point term and define the non-self-contracting (NSC) 
vertex : 

^"(x) = [e"(x)]^^^ + (all possible e contractions), (18) 

[C{x)f^^ = C{x) - (all possible ^ contractions). (19) 

The contracted vertices (second term of Eq. (fT8|) ) are treated as shifts of the lower order 
terms. These contractions, as they contain the tadpole diagrams, are typically UV diver- 
gent. We use the dimensional regularization and absorb it into the higher order LECs. In 
this way, the coupling renormalization can be done in advance, and one can substantially 
reduce the number of remaining one-loop diagrams for an arbitrary correlation function. 
Note that {[Cix)]^^^)^ = by definition. 
The two-point vertex is the easiest example: 

leixr'"' = e{x) - {eix))„ (20) 

which is applied to the 4-th term of Eq. ([7j), and in this case, the ^ contraction is treated as 
a shift of S in the second term of Eq. (I7j). Its UV divergence is absorbed into Lg. 

With the NSC vertices, Li terms in Eq. ( |T7I) . and measure term Eq. ( fT6i) together, we 
can express the low-energy effective action as 

Jd^xC = -^^Tt[M^Uo + U^oM] 

Nf 

+ J d'xWizfr ([d^M^r'^'i^) + K)' [6,e,.]^'^(x)) 

+5«(^o,0+5f (t^o) 

+<Sdiag + <S3pt + Sipt, (21) 
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where the first two terms are the LO contribution, and the perturbative interaction terms 
are given by 



S?\Uo.O = / d'x^Ti \(m\Uo - 1) + {Ul - 1)M) 



2F2 



1 TV 5(7 



■Ml 



Here we have used notations below, 



i-;^|E^(o,My2)-G(o,o,o) 



' Nf 



n p2 



X:(A(0, Mj) - A(0, Mfj/2)) - (G(0, M^, M^) - G(0, 0, 0)) 



7*-? — ^ 



1 - - 1 - 

A(0, Mj) + A(0, M]f)) + -A(0, M^, M 
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TV, 



-8 |L4E^?/ + ^5M, 



and (M|.)^ = + Nf/F^V with 



=14- 



0(0, Ml Mj.) - 8(L4 - 2L6) ^ - 8(^5 - 2Ls)Mf^ 

f 



In the last line of Eq. ( l2Tl) , we have 
1 / 



A(0,M2.) 



- ( |MJ A(0, Mj) - l&L,MlMl + ^) 



33\ 



^spt = J d X 
Sipt = I d^x 






372^3 

1 ^/ 



(M^Uo - UlM 



(22) 



(23) 



(24) 
(25) 



(26) 



(27) 



X 



(28) 
(29) 

,(30) 



but they do not contribute to the calculations in this paper where we only consider two-point 
functions of off-diagonal sources. We therefore simply ignore them in the following sections. 
We have also ignored trivial constant terms in the above expressions. 
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C. Pseudoscalar (and scalar) source term 



The pseudoscalar and scalar source terms are obtained by extending the mass matrix: 

M^Mj = M+tJix), (31) 

where the pseudoscalar and scalar parts are given by 

pix) = ^{J{x) + J\x)), (32) 
s{x) = '-{J{x)-J\x)\ (33) 

respectively. 

In order to keep a manifest and consistent counting rule, we treat j in the same way 
as the original mass matrix, i.e., 

J{x)^0{p''), J{x){U,-l)^0{p'), J{x){U, + Ul-2)^0{p'). (34) 

Note however that unlike the original mass matrix, ^/-derivative could isolate the matrix 
element of (f/o — 1), which could cause ambiguity in the counting rule of correlation functions. 
In fact, the leading contribution of the pseudoscalar two-point function is known to be (9(1) 
in the e expansion while it becomes one order higher, O(p^), in the p expansion. To avoid 
this problem, we consider every iTij-derivative multiplied by a factor ^irtimj 

(35) 

as a unit block of the calculation. This prescription keeps the counting order of the operand 
unchanged even after differentiation. Note that the unusual square root does not appear in 
the physical results since even numbers of derivatives are always required to give a non-zero 
correlation when i ^ j. The pseudoscalar two-point correlation, which is our target of this 
work, is then kept always at 0{p^) in an unambiguous way with arbitrary choice of the 
quark masses. 

Unlike the Lagrangian itself, we need to introduce an unphysical constant counterterm 
with a coefficient H2 |3|, 

- H, (^p)' Tr[(A^ + iJ)\M + (36) 
to cancel the divergence of the scalar operator at a finite valence quark mass. 
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SJ{x) 



Now let us collect terms linear in J and rewrite it in terms of NSC vertices at O(p^) 



ij \2 

m) 



X 



Nt 



4(2L8 + H2)Ml 
F2 



3F3 



2F2 



Tr[:rt(x)f/oe'(x) - e'(x)f/o^:r(x)] 



(37) 



where a term with the cubic NSC vertex [^'^]^'^*" is ignored since it never contributes to the 
two-point correlation functions. A new factor Z^p is defined by 



— 1 



E(A(0, Mj) + A(0, M]^)) + A(0, M^, Ml) 



Nf 



-8\L,Y.Mjf + L,M, 



(38) 



D. Renormalization 



In the above results, A{0,M^) and (5(0, Mf,M|) have the exactly same logarithmic 

divergences as the conventional p expansion since the absence of the zero mode do not affect 

the ultra-violet properties. In the same way as in js], we can thus evaluate their divergent 

parts by the dimensional regularization at D = 4 — 2e (taking e <^ 1): 

M2 /I 
IGtt^ ve 



A(0,M2) 
G{0,MlM^) 



„ . - + 1 -7 + ln47r 
1 IMl + M^ 1 



167r2 



N 



f 



^ + 1 - 7 + ln47r) +■■■, (39) 



where 7 = 0.57721 ■ ■ ■ denotes Euler's constant. As is the usual case, these divergences can 
be absorbed into the renormalization of Lj's and H2 as 



sub) 



li 

327r2 



^ + l-7 + ln47r-ln/i^„f,^ 



H2 = H^if^sub) - + 1 - 7 + ln47r - In/i^^^ ; , 



(40) 
(41) 
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where L^^figubYs and Hl^i^Hsub) denote the renormahzed low energy constants at the subtrac- 
tion scale ^sub and 



Nf ^ , ^ 

74 = -, 75 = — , 76 = o I o + I ' 77 = 0, 78 



1 

8' 



8 I 2 



7H2 



1 fN 



As a result, Scg, AZj^, Zp and Z]^ are kept finite, while Z^"' still diverges but it never 
appears in the physical observables. 

After this procedure, one can replace A(0, M^) by. 



A''(0,M2) 



16^ 



ln^+^i(M2 

f^sub 



(43) 



where gi denotes the finite volume contribution of which the zero-mode part is subtracted. 
It is well-known that there are two expressions for g: one valid for small ML ^ 1 ll|] and 
the other valid for ML ^ 1 , and their convergence around ML ~ 1 is discussed in detail 
in Ref. 3J]. Here we just note that on a L ~ 2 fm box, these two 



r \n^\<n^^ /T72 



MM' 



M^ 
IGtt' 



47r2| 



M'^V 



{\M\L > 2) 



(44) 



ln(MVi/2) - J2 



^2(n-l)^(„-2)/2 Q^i^ < 2) 



at n^'^^' 



1 (^-1)!' 

7 and n™"^ = 300 show a good convergence around the threshold \M\L = 2. Here 
Ki is the modified Bessel function and the summation is taken over the 4-vector a,,, = n^L^ 
with Lj = L (z = 1, 2, 3) and L4 = T. /3j's denote the shape coefficients defined in 



IV. ( CONTRACTIONS IN THE CORRELATOR 



We are now calculating a hybrid system of a matrix Uq and fields ^ whose partition 
function (with the source J') is given by 



dUo (det t/o)^ / exp 

U{Nf) JsU{Nf) 



J d^x (£ + 



(45) 



where we need to integrate over both fields. The integral over Uq, in particular, has to be 
non-perturbatively performed. Our strategy of this study is (i) to perturbatively calculate 
^ fields first, (ii) then to perform Uq group integrals. 
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Let us here define two notations 



dUo (detf/o)« ^^TrlMWo+U^,M] Q^^^^^ 

{Oi{Uo))uo = -r — ; , (46) 

'dUoidetUofe-^^^^^'^o^^'oM] 



^^'^^^^^ ^ j die^^ ''"^ Y.a^'i?[i^3{-al+M?,)i..f'^ ' ^^^^ 
with which any correlation function of Uq and ^ (we denote f{Uo,C,)) can be expressed as 

(/(^o,0) = ^ \ (48) 

where the interaction terms iS}*'''s are treated perturbatively. Noting Sj^\Uq,^) ~ 0{p) and 
iS}^'*(f/o) ~ C^(p^), the correlation function above at NLO can be divided into four parts: 

(/(f/o,o) = (/(f/o,or+(/(f/o,o)'°+(/(f/o,o)'°+(/(f/o,o)°s (49) 

where the superscripts 00,10,20,01 mean C(l), C(>Sj^^), C((5}^V) and C(5f ^), respectively. 
Namely, they are defined by 

(/(t/o,Or = ((/(f^o,0)c)^^, (50) 
(/(f/o,0)" - ((-'5l'^/(f^o,0>J^^^-((-5f^>J^/(/(f/o,0),>^„, (51) 



These notations are useful in the following calculation. 

In the rest of this section, we calculate the (■ ■ ■)g part using the Feynman rule, Eq. (ITOl) . 



A. Chiral condensate to NLO 

For a warming-up, let us first calculate the one-point scalar function {i.e. the chiral 



condensate) to the next-to-leading order 
diagonal matrix element of the source 



In this case, we consider a pure imaginary 



[J]ij = -i6iJjySyy{x). (54) 
14 



In this case, the source term in the Lagrangian is 
Cj = 



-y^s„,(a;)[f/o + Ul]^^ - i-j=^ 



[x)\U^i{x)-i{x)Ul\ 



E 



TO 



F2 



2F2 



s„„(x)[^of (x) + f (x)f/ot]f/^ + 0{p 



where the index v is not summed over. 

Now we can calculate the chiral condensate of the w-th valence quark as follows, 



m^{qvqv{x)) = 



nir. 



^^vv (^) 
/ Seflf 



Sw — 



^ 2 

Seff 



r2 \ \ 00 



{[Uo + uIu)u,-t.[az, 



E 

vv 



4(2L8 + H2)M^ 



2 ' 

lit; 



i?2 



(55) 



(56) 



where we have used {i{x))^ = and (l)c;o = (l)^ = 1- Note that {(jvQvix))^^ = {(jvQvix))'^^ = 
{qvqy{x)Y^ = to our order, which can be easily confirmed by a direct calculation using the 
fact {\Uq + f/d]t,t,)°^ = {\Uq + Ul\vv — 2Y^. The result is, of course, consistent with Ref. 34|. 



B. Pseudoscalar correlator 



Let us next consider the pseudoscalar source. In the calculation of meson correlators, 
we take a specific generator of the chiral group which has V1V2 and V2V1 (f 1 7^ V2) elements 
only. This choice corresponds to the charged pion or general kaon type correlators. Here Vi 
denotes the valence quark index whose mass is given by m^. . For simplicity, we omit "w" in 
the following: the indices Vi and V2 are denoted by 1 and 2, and their masses are expressed 
by mi and m2, respectively. Namely, we consider 

= ]j{.hi^j2 + 52i5ij)p{x), (57) 

where p{x) is a real classical number. 

The pseudoscalar source term in the Lagrangian then becomes 

Cj = ^(P-(x)+P-(a;)), (58) 
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where 



£f {[Uo\i2-[UlU) (l-AZ 



22 + -p^^l2 



'12 vl2,' v12n2 



^E([f^o]l^62(x)+6^(a:)[t/o^]^l 



2F2 



E[e'(^)]f''([^o]i.5,2-<5iat^o'].2 



(59) 
(60) 



Now we are ready to calculate the pseudoscalar-pseudoscalar (PP) correlator, 



mim2{P{x)P{0)) = 2mim2 



6 



6 



mim2 



Z(0) 5'p{x) (5p(0) 
1 



Z{J) 



1 



p(x),p(0)=0 



2(P^^(x)P^i(0)) + -{P'\x)P'\Q)) + (1^2) 



(61) 



where an overall factor of 2 is introduced to compare with the corresponding lattice connected 
diagram. Note that the procedure Eq. ( 135|) is performed but the factor 17111712 will be omitted 
for simplicity in the following calculation. 

Although the number of diagrams we need to calculate is substantially reduced by us- 
ing the NSC vertices, our calculation is still tedious because of the off-diagonal elements of 
Uq in the source term Eq. (|59l) . which produces various unusual channels in the correlator. 
Every step of calculation is, however, rather straightforward as in the conventional p ex- 
pansion, except for the use of the A^(f/o — 1) ~ 0{p^) rule. We therefore skip the details 
of the calculation in the main text here. Instead, we summarize several useful formulas 
for the computation in Appendix \K\ and present each piece of {P{x)P{0))^^, {P{x)P{0))^^, 
(P(x)P(0))2° and (P(x)P(0))°^ in Appendix [Bl We also use the technique in Appendix [Dl 

After relevant one-loop integrals over ^, the pseudoscalar correlator is given by 



(P(x)P(O)) = (P(a;)P(0))°° + (P(x)P(0))i° + (P(x)P(0))2° + (P(x)P(O)) 

4 /il + /i2 



01 



^cff 



Ob 



+Y(A^fi-AZi)C°^+^ 



(Z-(Z-)2)2C^A(x,M-) 
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+Cf, [A{x, Mi\) - A(x, M^,)) + Cl (A(x, MI2) - A(x, M\ 
+ E Ct, (A(x, MJ.) - A(x, M',,)) + 5: 4 (A(x, M^^J - A(x, M 



+C'>G{x,Ml,Ml,) 

+Cl^ (G{x, Ml^, Ml^) - G{x, M\ 



+Cl, (G{x, Ml, Ml^) - G{x, Ml,, M 



22, 



(62) 



where 



Uo 



706 



+ 



Uo 



(63) 
(64) 

(65) 



Nf 



Nf 



C ^ (i[Uo]n + [mh2)i[Uo]22 + [m]u) + Y.[Uo]iAUl]n + Y.[UoUU'^^^^ 

\ j^l i^2 



Uo 



^ ( 2([7^]n + [nu - y: - E f^]"''^]''^ 



C., = 2(([f^oy +([f/oy )c/o + ^^^^ + 2(m, - m,)2 ^ 



(66) 
(67) 



= \FoJii['^oJji/f/o H ^ 



mj — mj {rrij — rriiY 

C ^ - (([f/0]l2 + [f/J]2l)([f/o]21 + [f/J]l2) + ^([f/o]l2 + [UlUf + ^([f/o] 



(68) 
(69) 



21 T~ [<-^oJl2 



ci ^ \m,h + [ul].,f)uo + 



{{\n^ + mMUoh + [ulV,))uo 



(2[7^].,[7^],, + ([7^].,)^ + ([7^y^)^„ 

2(mj — mjY 



where we have used a notation 

7^ = 7Wt(t/o - 1) + (f/d - 1)A^. 



(70) 



(71) 



(72) 



One should note that many unusual channels appear in Eq. ( 162|) . which is a quite un- 
natural situation when just a single particle propagator is expected. However, one will find 
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in the next section, many of them actually disappear, or many of the coefficients C's vanish 
after integration over Uo ^. 



V. ZERO-MODE INTEGRALS 



The zero-mode's contribution to the so-called graded partition function of n bosons and 



m fermions is analytically known 



47 



(( IN ^ det[/i] jQj^j_i{lii)]ij=i^...n+m /„„x 

""--^^^^'^^ n ■>.=i(/i| - i{-tr=n+M - /^D ' ^ ^ 

in a fixed topological sector of Q where /ij = rrii'EV. Here J^'s are defined as j7Q+j_i(/ij) = 
{-iy~^KQ+j^i{fii) for i = 1, ■ ■ - n and Jq+j^i^fii) = Iq+j^iifii) ior i = n + 1, ■ ■ ■ n + m, 
where Ii, and Ki, denote the modified Bessel functions. Partial quenching is completed by 
taking the boson masses to those of the valence fermions at the very end of calculation. 

Exact group integrals of various matrix elements over Uq can be calculated by differenti- 
ating the above partition function. The most basic pieces are 

1 d 

= 7;([Uo]vv+[UXv),. = 1™ — InZi^i^^ (/ib,/i„,{/i,ea}), (74) 

1 

4 



Q2 

lim r2:i^i+^^,(/ife,/i^„ {/isea}), (75) 



d d 



where /if,^ denotes the bosonic spinor mass and {fisea} indicates a set of sea quark masses 
(normalized by T,V). Note that V^-^y^ differ even when m^,^ = = 

In Ref. 16|, more non-trivial matrix elements are calculated in terms of the above iS's and 
P's using the left and right invariance of the group integrals. Their results are summarized 
in Appendix [Cl 



^ The readers might wonder if the integration over first is then inefficient. But if we perform Uq integrals 
first, we need much more tedious computation over Uq than what we wih see in Sec. |Vl which do not 
disappear until ^ integration is completed. We thus believe our order of calculation is easier. 
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Now we can simplify C"s in terms of iS's and P's. Note here that for the leading con- 
tribution, namely for and C^, we need to use Sefr instead of S in the arguments. We 
distinguish them by putting a superscript "eff" like /if^(= mjSeffV^) and Sf^ . The results 
are summarized below. 

4 



/;Cff _|_ 



0, 

2 + 

2mim2 
0, 



(5f + 5f ) 



cff 
12 



eff cff / ' 



2mi(5i - 1) + 2m2(52 - 1) + E 



mi — TTij 



E 



m2 — mj 



4 



{Si - S2) . 



(77) 
(78) 
(79) 

(80) 

, (81) 
(82) 
(83) 
(84) 



/il - /i2 

Note that we have used mi{Si — 1) ~ 0{p^) 

Since the term contributes only in the p regime, one can substitute the perturbative 
expression to Si 



42 



Nf 



Si 



1-E 



+ 



+ 



(85) 



and obtain 



mim2C^ = 4mim2 



(mi + 1712) 



(86) 



2/ii/i2 

Noting mim2C^ = 4mim2 + 0{p^), the 5th term of Eq. ( l62l) can be absorbed into the 4th 
term (namely, term) by shifting the meson mass as 



N 



'12 



(87) 



F2\/ ' 2/ii/i2^' 

We recall that an unexpected term is found in the definition of but it is now 
canceled out. 

Thus the result can be expressed in a simpler form. 



cff 



(P(x)P(0))q = T:\Zff {Z''^' 



M) 



Sf + SI 
AH + /i2 
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where 



+ 



712^2,' ^12-14 



l + V 



cff 
12 



eff efr 
PI P2 , 



F"^ /il - /^2 



(88) 



(89) 



VI. RESULTS 



A. Pseudoscalar correlator at fixed topology and in = vacuum 



Let us take the zero-mode projection, or integrate Eq. flHHjl over three-dimensional space 
(See eq.dXT])). 



VV{t,mi,m2) 



d'x{P{x)P{0))Q 



12\2 o I ^ ^ ^12 



2 

7l2 



^\ + At2 



5f + Sf 



cosh(Mg(t-T/2)) 
MVf y Mi^2 sinh(Mg^/2) 



■ off, cff 
PI P2 , 



2/ilAi2 



2S^ S\ — 



111 - /i2 

which is more useful to compare with lattice QCD results, where 

r,,(t) = / £xG{x,MlMl). 



(90) 



(91) 



This is our main result in this paper valid for an arbitrary number of non-degenerate flavors. 
It is also important to consider the correlator in the ^ = vacuum. 



PP(t,mi,m2)e=o = / c/'x(P(x)P(0))e=o 



T?{ZlZ^l f 



F^{ZfY 2 



(^^?f )^=o + 



(QVo\ cosh(Mf=°(t-T/2)) 
liflif I Mf2=°sinh(Mf2=°T/2) 



1^1 + 1^2 

2Y? (iSi)6i=o - 



(iS2)6l=0 



1 + (^12 )e=0 + ^ 'eff'^cff 
V PI P2 



61=0 



1 + 



(gVo^ 

2/ii/i2 ^ 



F2 



/^l - /i2 



-ri2(t), 



(92) 
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where {M^f^ ^ M^. {z% + 



The summation over topology, 



EQ0(Q)^o%({/^f}) 



(93) 



can be, at least, numerically performed using the analytic expression for which 
is finite. For small Nf cases, simple analytic forms are also known 50| • Note in the p regime, 
that we can easily calculate {Q^)e=o = = 'mllV = XtV where m = 1/I]/(V'^/) ^^"^ Xt 
denotes the topological susceptibility ^. 

As seen above, we find a constant contribution in the pseudoscalar correlator in addition 
to the conventional cosh function of time t. This constant term is indispensable for keeping 
the result IR finite and giving a smooth interpolation between the e and p regime limits. 

B. Check in the p regime and e regime limits 



Let us confirm whether our above formulas recover the conventional p expansion results 
when both of mi, ar e large (or mi,m2 ^ l/SV). In that limit, we can use (see 



Appendix O and Refs. [421, ki ) 
1 



5,. 



Oip'), 

1-E 



V 



1 , Q' 

1 



12 



S1S2 ~ 1 - X] 



0{p' 



E 



+ 0{p' 



{Q^)e=o ~ /i- 



(94) 
(95) 

(96) 
(97) 



Here one should remember that in the conventional p expansion, Z factors are expressed 
not by A(0,M2) but by A(0, M^) = A(0, M^) + 1/MV. To take this into account, it is 
useful to redefine the Z factors, 



r/Yl 



7\2 



1 



4/ii/i2 
1 



1 / 1 



1 



(98) 
(99) 



In the p regime, the LO calculation of xt is enough in this work. See Refs. 



51 



52 



for the NLO correction. 
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The well-known result in the p expansion is then precisely recovered, 

VV{t,m,,m,)e=o - ^2([^^2]^)2 Mfr° sinh(Mfr°T/2) ' 
where M^^° = Mi2[Zl]]p. Note that the constant term and ri2(t) term rapidly vanish as 
mi or 7712 grows. We also confirm that our result at fixed topology agrees with the one in 



the p expansion 



Next let us consider the e regime limit, where both of the valence masses are near the 
chiral limit, mi ~ m2 ~ 1/TiV. In this case, one can expand the hyperbolic cosine term in 
the meson mass as 

cosh(M(t - T/2)) _ 2 



Msinh(MT/2) M^T 



2Thi{t/T) + 0{M'), (101) 



where 



and obtains 



Ceff I ccff OV2 C C 

+L'^l.%^ - %-^ri2{tl (103) 
/if + ixf Ui - /i2 

which is consistent with the result in the e expansion (Ref. 18|). Note that we have used 
C. When m2 is large 

One of our main interests in this work is to consider when one valence quark is always 
large, in the p regime: m2T,V ~ (9(l/p^). Namely, we consider the chiral limit of the 
kaon-type correlators in a finite box. 

In this case, we can perturbatively treat (see Appendix [C]) 

^ - ~ Oip'), (104) 



Hi + 1X2 fJ'2 

S2 + ^ + 0{p'), (105) 
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and the correlator in that hmit is 
VV{t,mi,m2)Q 



cff 



Y /i2 + /i/ 



X 



cosh(Mi^(t-r/2)) 



(107) 



The result in the 6^ = case is obtained by replacing Q with (Q )6i=o and 5^ with (5^ )e=o- 
One can see that the overall factor (and therefore the calculation of the decay constant 
Fk) still has a large finite volume correction from the zero- mode integration, while the meson 
mass [Myi here) has a rather small perturbative correction. 



D. Origin of the G[x, Ml^.M^^ term 

The third term in Eq. fl90|) becomes significant only when both of mi and m2 are in the 
e regime. Here we consider the origin of that term. 

Although non-perturbative integration of the zero-mode is supposed to be the most re- 
liable way of calculating the finite size effects near the chiral limit, it obscures the physical 
meaning as propagation of the pions. Let us here go back to a perturbative picture in the def- 
inition of Eq. fl70|l and express the corresponding correlation function using Appendix [D] and 
putting labels "(x)" and "(y)" to explicitly show where the original operators are located. 
For example, the first term of Eq. ( ITOl) is expressed by 

(([f/o(a:)]i2 + [f/o^(x)]2i) ([t/o(?/)]2i + Vlm^^))^^ 

= - Jj i{[U^)]i2[Uy)hi)u. + {[U^)hi[Uy)]i2)u.) 

(([eo'(^)]i2 + mx)himo{y)hi + my)]i2))^^^ 

+ J^i ((([eo(x)]i2 - [Ux)hi){my)hi - my)]i2))uo 

+mix)]i2 - mx)himo{y)h - uy)]i2))uo) + ■ ■ ■ .(los) 

With this perturbative picture of the zero-mode, the term can be expressed as 

c'G{x - y, Ml,, Ml,) = -1 {m{x)h + mx)himo{y)]2i + my)]i2))^^^ 

X mx)]u[ay)]22 + Ux)]22m]u))^ + ■■■■ (109) 

It is then obvious that this term is originally a three-pion-state propagator which is sup- 
pressed in the ordinary p regime. As the system enters the e regime, however, two of their 
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zero-mode's contributions are non-perturbatively enhanced and it becomes an NLO contri- 
bution. 



VII. USEFUL EXAMPLES 



In this section we present two specific examples in the A'^^ = 2 (with degenerate up and 
down quarks ) and 2 + 1 (with up, down and strange quarks) theories, which are useful to 
analyze lattice QCD results simulated in finite volumes. In the formulas below, we denote 
the sea quark masses by m„ = ma = m and (/i = mHV and fis = rrisTiV). 

We consider two-types of the pseudoscalar correlators: the pion-type correlator whose two 
valence masses are degenerate, mi = m2 = ruy {fi^ = m^'EV), and the kaon-type correlator 
for which we take m2 always to be in the p regime (see the general formula Eq. fll07p ). 



A. Simplified G(x,M2,Af2) 

For small Nf, we can simplify the G{x, Mf^, M22) (or ri2{t)) term. Since it con- 
tributes only when both of mi and m2 are in the e regime, it is sufficient to con- 
sider the pion-type correlator case with mi = m2 = rriy. The result was already pre- 



sented in Ref. 



5l[ |. except for the presence of the zero-mode part : G{x, Mf^, M,^ 



2 1 
22) 



G'(x,Mfi,M|2) + l/{VMf^M^2iT,f'^/M'ff)), which does not affect the coefficient of each 
term. Here we just present the results for the A'^^ = 2 and 2+1 cases, 

^ 'A(x, Ml) + (Ml - Ml)dM^^A{x, Ml)] {Nj = 2), 



G{x,Ml,M[ 



vv) 



1 {Ml 



2 {Ml - M2)2 
{Ml - Ml) {3(M, 



A(x,M, 



1 {Ml 



2 

vv 



Ml) 



V 2{Ml-Ml)^) 
{Ml -Ml)} 



-^\Hx,Ml 



2{Ml - Ml) 



9AeA(x,M2j 



{Nf 



;iio) 



where Ml = 2m^S/F2, Ml = 2mS/F2, and = (2m + 4m,)S/3F2. Noting that 
A(x, Ml) rapidly converges to —1/MlV for large |x| and remembering that the correspond- 
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ing term contributes only when ~ O^p"^), it is sufficient to consider (See also Eq.f lAip .^ 



cosh(MQ(t-T/2)) _ 1 
2M£ sinh(M£T/2) {M^vfT 

cosli(M£(t-T/2)) 1 



6M2T 



{Nf = 2), 



(111) 



.2M£sinli(M£T/2) {M^vfT _ 
Here we have used an additional assumption that the valence pion mass is not taken very 
differently from the physical pion mass and the 0{Ml^ - Ml) contribution is ignored. The 
only exceptional case : ^ M^^ will be discussed later. Note that we have replaced the 
tree- level mass by the NLO mass for later convenience (the difference is NNLO.). 



B. Nf = 2 and 2 + 1 flavor results 



Using Eq. f illip , the pion- type correlator can be expressed in a compact form, 



where the valence pion mass is given by 



cosh(MQ(t-T/2)) ^ 
m£ sinh(M£T/2) 



m2. 



M^^ZZ 1 



and 



pp 



S2 (ZX5'Z]jr)4 1 
F2 (Z^^)2 2 

25: 



cff 



95: 



cfT^ 



D 



-3 ^cff 
off 



PP 



3 ^cff 





cff 



i + p: 



cff 



cff ; ' 
dSf^ 



25: 



cff 



1 + v: 



cff 



ds: 



cff ^ 



^955\ 



cff 



Here we have used E^s = ^{Z]^Z]?f + O(m^) and 



lim 



5i — ^2 



(112) 



;ii3) 



;ii4) 



{Nf = 2), 



(iV; 



(115) 

2 + 1), 

(116) 



mi->m2=m„ _ 

It is also possible to simplify the kaon-type correlator (here we choose the second valence 
mass to be the physical strange quark mass: m2 = rris in the 2+1-flavor theory) as 



r(, ^ - -v-vf, ^ i^Q cosh(MQ(t-T/2)) 



;ii7) 
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where the valence kaon mass is given by 



Ml 



M^sZl^ ( 1 + 



(118) 



and 



^pp - p2 ^zff 2\:^^^ y 



L+^^ 



(119) 



The result in the 6 = Q vacuum is obtained by simply replacing with ((5^)e=0; 
with (5^''^)6)=o and with (P^![^)e=o in the above formulas. 

Using a notation for the renormalized logarithmic term which is given in Eq. (143|) . the 



explicit forms of Z factors 



53 



54| . Seff/S, and "D^^ (see Appendix are given by 



Nf = 2 case : 



M 



ryvv 



^off 



1 + 



2F2 



^A'^(0, Ml) + i(M^,, - Ml)dM^^K^{Q. Ml 



-IQ{LI - 2Ll)Ml - %{Ll - 2Ll)M; 



2F2 
1 



2A'(0, (AC + Mi)/2) - i(2L\Ml + L^M, 



2 ^ 

VVJ 



(120) 
(121) 



2A^(0,M2/2)-- 



/3i 



M2 



167r2 



lnyi/VL-/32 



32L^M2 



(122) 



det 



1 



-9^„(/i^A"Q+l(/i^)) /i^/Q + l(/i^) /ijQ + l(/i) jQ(/i) 

dt,,{fJ^lKQ+2{fiv)) flllQ+2{f^v) /i2/Q+2(/i) /i/Q+l(^) 
-^M.(/^^^Q+3(/i^.)) fillQ+3M H^Iq+M /^^^+2(/i) 



det 



/Q(/i) /X ^/Q-l(Ai) 



(123) 
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det 



-df,,ifivKQ+i{fi^)) (9^„(/i^/Q+i(/i^)) /i/Q+i(/i) Iqifi) 



( 



det 



/Q(/i) /i ^/Q_l(/i) 



(124) 



1 case : 



M 



1 + 



2F2 



1 (M^ - M2 



2\2 



+ 77 



6 - M2)2 
1 (M^^ - Ml 



1A^1(M|_ 
3 V 2 



- M^) - {Mi - M^)}aM.„ A^(0, M, 



6 - M2) 
8(L^-2L^)(M2 + 2M,^)-8(L 



vvJ 



ryvs 



1 + 



2F2 



I Ml 



^Ml 



M2 _ 



M2 



I Ml 



2L^)M, 



M2 - 



ray 



-8(11 - 2Ll){M^ + 2M^) - 8(1; - 2L^)M, 



(125) 



(126) 



ryvs 



1 - 



^efr 



2F2 



2F2 



2A^(0, {Ml + Ml) 12) + A^(0, Ml) - 8{L\{Ml + 2M|.) + L^^M; 



A^(0, (M„2„ + Ml)l2) + -A'-(0, M,2J + A^^O, Ml) 



|-2 N 



+ 



1 / 1 A/2 

n ^ 

3 \ 2MI, 
1 I Mi - Ml 



M^V 

1^) ^'("-^'ii 



3 ML 



1 (M, 



2 

vv 



M2 

M2) 



1 

2 



1 
4 



> ^ A'^(0,M2 



M2 



M2 



{ml - Ml 



12 {Ml - M2) 
-8(LKM2 + 2M^) + rgM, 



-{Ml-Ml}dMsA'{0,M, 



' 1 



2A'^(0,M2/2)+A'^(0,M,y2) 



(127) 



(128) 
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1 

3 



Ml{Ml 



2M4 



{Ml 



Pi 



2M4 



2M2 



In V^/^L - /32 ) - 16L^(M2 + 2Ml] 



(129) 



5„ 



(^2 _ - /i2; 



det 



-9^„(/i„i\:Q+l(/i^)) /i„/Q+l(/i^,) /Q(/i) /ijQ+l(/is) 

dt,,{filKQ+2{fiv)) fillQ+2M /i^/Q+2(/i) /i^/Q+2(/is) 

-d/^Af^'l^Q+sM) Ullq+sif^v) /i2/Q+3(/i) /i2/Q+2(/i) /i^/Q+3(/is) 

V '9^.(/^^^Q+4(/^»')) /^^^Q+4(/it.) /^%+4(/^) /i^/Q+3(^) ^1^+4 (/is ) y 



det 



^ /Q(/i) /i^^/Q_l(/i) /Q(/i,) ^ 



,(130) 



X 



(/.2-;,2)2(^2_^2) 

-9^„(/i^,KQ+i(/i^,)) df,^{lxjQ+i{fi^)) filg+iifi) Igifi) fislQ+iifis) 

^^^MKQ+2M) df,^{fillQ+2{f^^)) /i^/Q+2(^) ^^Iq+M ^s^+2(/is) 

-dt,,{f^lKQ+3{fiy)) df^^ifillq+sM) fi^Iq+sif^) fi^lQ+2{f^) fillq+sif^s) 

\ '9^„(/i^^Q+4(/i^)) df,^{^tlQ+4{fi^)) fi^lQ+iifl) fi^lQ+3{f^) /i^jQ+4(/is) y 



det 



det 



^Q(/i) /i ^^-l(/i) ^(/is) 



+ 



+ 



2Ai„ 



Here we have used explicit expressions for G(0, M^, M|)'s shown in Ref. [511 and 

/3i 



hm A70,M2) 

A/-5.0 

hm aM2A70,M2) 

M-5-0 



1 



IGvr^ 



(131) 

(132) 
(133) 



For and 1^^,^, at degenerate up and down quark masses, we have used an expansion 
(/i + A/i)"/Q,(/i + A/i) = fi^Iaifi) + /iA/i[/i"^-'^/ct„i(/i)] + C((A/i)2) for any a, and a similar 
expansion for K^^s. Note that 5*^^ and are obtained by simply replacing S with Sgfi in 
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the above formulas. 



C. When M„„ < 



In Eq. fillip , we have neglected a term proportional to M^^, — M^. One might, however, 
encounter the case where one wants to reduce the valence quark mass to the very vicinity 
of the chiral limit while keeping the physical pion mass at the p regime. In such a case, a 
partial quenching artifact is enhanced as a double-pole contribution and one has to add the 
following contributions to the pion correlator, 

"cosh(M^^(t-T/2)) 1 



A7r(t, 



-FppdM^^ 



where 



vv 



pp 



vv 



Ml 



Mi 



Ml 



3 Ml, 



M2 



2M,, sinh(M,,T/2) M^T ^ 



(Nf = 2), 



(Nj = 2+1). 



(134) 



(135) 



D. Masses and decay constants 



In this subsection we demonstrate how to extract the masses and decay constants of 
the pions (and kaons) from lattice QCD data using our formula. We plot in Fig. [T]the pion 
correlator Eq. flll2p (normalized by S) at several different quark masses. We take rriud = 'n^'v 
in all cases. In the plot, the strange quark mass is fixed at rris = 111 MeV, and the topological 
charge is fixed aX Q = 0. We choose the finite box size as = L^T = (1.8 fm)^ x (5.4 fm) and 
the boundary condition is periodic in all directions. For the inputs, we use one of the latest 
lattice QCD results for the chiral condensate and the pion decay constant, E = [234 MeV]^ 



(in the MS scheme at 2 GeV) and F = 71 MeV from Ref. 41[. For the other low energy 
constants, phenomenological estimates from Ref. [3j, Ll{770 MeV) = 0.0, Ll{770 MeV) = 
2.2 X 10-3, L^(770 MeV) = 0.0, and L^(770 MeV) = 1.1 x lO'^ are used. 

As the first step of the analysis, one should identify the presence (or absence) of the 
constant term Dpp, which is a signal of entering (or leaving) the e regime. As shown in 
Fig. [21 it is a rapidly decreasing function of the quark mass. Since this constant comes from 
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the zero-mode part, it is essentially controlled by the chiral condensate. Using lattice QCD 
data for S (or Scfr) or taking time derivative of the correlator, -Dpp can be subtracted. 

Next, from the remaining cosh function part, the meson masses can be determined. In 
Fig. 131 we plot the quark mass dependence of the pion mass squared divided by the quark 
mass: {M^=°/^=°)^ / {2mud) and that for the kaon mass: (M^=°/^=°)V("^«d + ^ri^). Here the 
same inputs shown above are used. The ^ = results here and in the following are calculated 
via Eq. fl93l) truncating the sum at \Q\ = 20, which already shows a good convergence. For 
the pion mass, 10-20% deviation from the infinite V result (thick curves) is found near the 
chiral limit while the kaon mass suggests only ~1% finite volume effects. Note that there is 
no contribution from the zero-mode to the meson masses at Q = (See Eq. fl89l) .). 

Finally let us discuss how to determine the pion decay constant from the coefficient Cpp. 



It is not difficult to check that a naive conventional definition = y AmlCpp/ {MmY 
its counterpart in the 9 = vacuum F^^^ = {F^)g=Q actually leads to the right infinite 
volume limit = ^^^""1^=00 as V increases. It is also the case for the kaon decay constant: 



= \J {rriy + nisYE^p/ {M^sY (or F^'^) converges to the infinite volume limit of Fk = 
^^^■^1^=00- Note however that the curves in Fig. H] show a considerable deviation (~ 50%) 
as the quark mass is reduced, which is a typical consequence of the non-perturbative zero- 
mode integrals. Unlike the meson masses, not only the pion decay constant but also the 
kaon decay constant receives a large contribution from the zero-mode. These zero-mode 
integrals are again controlled by the chiral condensate, and therefore one should in principle 
be able to subtract this part using lattice QCD data for S (or Scfr). Once the zero-mode 
part, - 1 + tJ^w - S or 5f (l - - ^ + g) + i - 1, is subtracted, one 



obtains F^ = FZ^p or = FZp^, which have a much milder volume dependence (at most 
a few % level) as shown by the dotted curves in Fig. |H 

We emphasize that the accuracy of our calculation is NLO even though the zero-mode 
contribution is partly treated to all-order. It is interesting to compare our results with the 
conventional finite volume formulas in the p expansion since higher order loop calculations 



55j for the latter. In Figs. E] and El we plot our results 



are available a la Liischer formula 
for 

M^/^=°(L) _ F?/^=°(L) 

"^^^^^ = M.ML = oo) - = F.ML = oo) - ^'^'^ 

comparing with those in the two-loop (and one-loop) calculations in the p expansion by 
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a. 



Colangelo et al. [7|]. The same inputs for S, F and Lj's above are used. For the other higher 
order LECs, the values given in 7| are used. 

Our formula at one-loop (denoted by i exp.) in the 9 = vacuum is drawn by the 
solid (T =5.4 fm) and thick (T =7.2 fm) curves while the dotted curves (T =5.4 fm) show 
the results from which the zero-mode contribution is subtracted. Note that even in the 
region M-j^L < 2, our formulas are finite while the p expansion (dashed curves) results 
show an unphysical divergence. For M^L > 2, on the other hand, we observe that our 
result is consistent with the p expansion. It is, in particular, remarkable that our one-loop 
result is closer to the two-loop formula rather than one-loop in the p expansion. In order 
to understand whether this is a just coincidence or can be explained by the effect of the 
zero- mode resummation, a further study in the limit of T — oo, which enters another 
regime (the 5 regime |56l-l62|). is needed. 



We have observed that, as the quark masses decrease, the pseudoscalar correlator in a 
finite volume is largely distorted from the form in the infinite volume limit because of the 
zero-momentum mode fluctuation. By a careful removal of its contribution using the ChPT 
formulas, however, we can obtain a milder volume dependence, which makes it possible to 
extract the V oc limit of the meson masses or decay constants. 
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FIG. 1: The Nf = 2 + 1 ChPT prediction for the pion correlator 7r(t,m^)Q (normahzed by S) at 
Q = and = rriud- The finite periodic box size is V = L^T = (1.8 fm)'^ x (5.4 fm). We use 
rus = 111 MeV, S^^(2 GeV) = [234 MeV]^, F = 11 MeV, ^^770 MeV) = 0.0, L^(770 MeV) = 
2.2 X 10^3^ L^(770 MeV) = 0.0 and Lg(770 MeV) = 1.1 x 10"^ as the inputs. 




20 40 60 80 

(MeV) 



FIG. 2: The quark mass mud{= t^v) dependence of —Dpp (normahzed by E) at Q = 0. The same 
inputs as Fig. [1] are used. 
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FIG. 3: The up-down quark mass mud{= f^v) dependence of (M, 



Q=o/e=o^2 



Y/{2mud) (upper panel) 



and {M^s '^)'^ /{rriud + rris) (lower) is plotted at different volume sizes. The same inputs as 

Fig. [1] are used. 
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FIG. 4: rnud{= it^v) dependence of the pion (upper panel) and kaon (lower) decay constants 
pQ—O/d—O pQ—o/d—o different volume sizes. The curves with the index "zero-mode subt." 
denote or F^. See the text for the notation. The same inputs as Fig. [T]are used. 
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FIG. 5: Comparison with the p expansion results a la Liischer formula 



551]. Our new ChPT 



calculation {i exp.) and the p expansion {p exp.) result from Ref. |7l] for Rm^, (top) and Rmk 
(bottom) are drawn (note that one-loop correction in the p expansion on Rmk zero). The same 
inputs as Fig. [T] and those given in Ref. {t] for the other higher LECs are used. The M^-L = 2 
(thick) curve is drawn using the two- loop result. 
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FIG. 6: Comparison with the p expansion results [Tj for —Rf^ (top) and —Rfk (bottom). 
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VIII. A SHORT-CUT PRESCRIPTION 



We have performed a complete calculation to obtain the general form of the pseudoscalar 
correlation function in Eq. flUU]) . which contains a conventional cosh function as well as a 
constant term and a contribution from 3-particle states. 

It is no surprise that the constant term appears since the correlator in the conventional 
p regime shows an unphysical infra-red divergence in the chiral limit. To remove this diver- 
gence, the zero-mode or the constant mode contribution is indispensable. 

With this observation we find that the result in Eq. fl90l) is obtained by an easier pre- 
scription below. Starting from the conventional p expansion formula in Eq. fllOOp . 



1. Replace the Z factors with those from which the zero-mode contribution is subtracted, 



namely, [Zp]p and [Z]^]p with Zp and Z 



M- 



2. Replace 



cosh(Mf=0(t - r/2)) 
M^=o sinh(M^=0T/2) 



with 



cosh(Mg(t-T/2)) 2 
mJ sinh(M;jT/2) ~ (Mg)2T' 



(137) 



3. Multiply a factor coming from the exact zero- mode integrals, which can be read off 
from the coefficient of the t dependent term or the 2Thi{t/T) term in the e expansion 
resuh. In the case of Eq. ([90]), it is ^ ( 1 + Pjf + ^^J^ ) obtained from Eq. f lT03l) . 



2 i -12 

Note that the NLO condensate Seff, which contains chiral- log terms, should be used 
instead of the bare value S. 



4. Add the constant and ri2{t) terms if they exist in the e expansion. 

In fact, in a similar prescription, it is not difficult to obtain (a conjecture for) the 
axialvector-pseudoscalar and axialvector-axialvector correlators: 



AV{t,mi,m2)Q 



d'x{Ao{x)P{0))c 



12 ryl2\2 



V 



off 
12 



cff cfT 
PI P2 , 



1 + 



-j:{zyz 



M, 



-cff 



,cff 
12 



2yUi/i2 ) 



-Ml^^^^^-{dMir,2{t) + 9,,/|^ri2(t)), 

/il - /i2 " 22 



sinh(Mi^2(i-^/2)) 



sinh(Mi^2^/2) 



1 + 



2/il/i2, 



t 

T 



(138) 
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AA{t,mi,m2)Q 



d^x{Ao{x)AoiO))Q 

fni + 7712) (Si 



sf 



cosh(Mg(t-T/2)) 
sinh(Mi^T/2) 



12\2 



T 



^5? 



cflF 



+ 



T 
2^ 



off 
12 



cctt 



2/^1/^2 



1 + 1?' 



,cfr 
12 



, cff, eff 
PI P2 , 



Nf 



, cff, eflf 
PI P2 , 



where 



-1 



1 



(139) 



(140) 



^=(,,,,,,3) 4sinh2(y'|q|2 + M2 T/2) M^T^' 

which is UV finite (and of course IR finite as well) and can be thus numerically evaluated. 
We confirm that Eqs. (11381) and Eq. (I139p indeed converge to those in the p expansion 51 1 



16 



near the chiral limit. The above 



for the larger masses and those in the e expansion 
prescription thus achieves at least a smooth interpolation between the e and p regimes. Note 
that the e regime result is not found in the literature for the AV{t, mi,m2)Q correlator. We 
present in Appendix [E] our own calculation. 

Furthermore, we find a more non-trivial evidence that supports our prescription: the 
axial Ward-Takahashi identities 
d 



dt 



AV{t, mi, 1122) Q = {mi + m2)VV{t,mi,m2)Q 



12 V12N2 



M) 



eff eff 
PI P2 , 



1 + 



T 

Si — S2 

fll - fl2 



{sf 



Ccff 
'^2 



2fllfl2 J 



Mgcosh(Mg(t-T/2)) 
sinh(Mi^2^/2) 



2/il/i2, 



-2E 



and 



d 



- —AA{t,mi,m2] 
at 



M,Vi2(t) + 0(/), 



(mi + m2)AV{t, mi, m2)Q 



:i4ii 



T.{ZYZ]jY{mi + m2) 



Sf + Sf sinh(Mi^2(^ - T/2)) 
2 sinh(Mi^2^/2) 



+ 0{p^ 



(142) 



are precisely satisfied. Here we have used 



{Sf 



ccff 
'^2 



l+V 



cff 
12 



, cff , cff 
PI P2 , 



2/il/i2, 



t 

T 
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,cff 



1 



2/ii/i2 



sinh(M{^(t-T/2)) 



(143) 



2 sinh(Mi^T/2) 
it, and 

z = l,2). (144) 



which is vahd up to a higher order contribution near the chiral hmit, and 

Our resuhs in Eqs. (!90|) . fll38p and fll39p not only smoothly connect the e and p regimes but 
also keep the symmetry of the theory even in the intermediate region. 



IX. CONCLUSION 



With the new perturbative scheme of ChPT proposed in Ref. 3J], we have calculated 
the two-point correlation function in the pseudoscalar channel. The counting rule for the 
computation is essentially the same as in the conventional p expansion (except for the 
additional rule for the mixing term of the zero and non-zero modes) while some of the 
zero-mode integrals are performed non-perturbatively as in the e expansion. 

As seen in Eqs. flQOl) and Eg. (11 121) . the correlator is expressed by a hyperbolic cosine 
function of time t plus an additional constant term as well as a non-trivial contribution from 
three-particle states, which smoothly interpolates the p regime results and those in the e 
regime. 

The presence of the constant term in the correlator was known as a remarkable feature 
of the e expansion. We have found that this constant plays an essential role in canceling the 
unphysical divergence coming from the cosh term in the p expansion and keep the correlator 
always IR finite. 

Giving examples for the Nj = 2 and 2+1 theories, we have proposed a new method of 
determining the meson masses and decay constants from lattice QCD data obtained in a 
finite volume. Once one has a good control of the chiral condensate E, and therefore, of 
the non-trivial coefficients Cpp, Dpp and Epp in the correlators Eqs. (11121) and (11171) . the 
zero-mode contributions can be subtracted and the remaining meson masses (see Fig. [3]) and 
decay constants (Fig. Hj) show a much milder volume dependence. Our results will be useful 
to precisely estimate the finite volume effects in lattice QCD data for the pion mass and 
kaon mass Mk, as well as their decay constants and Fx- 

From our calculation we have found a short-cut prescription as shown in Section IVIIIi 
According to this greatly simplified scheme, we have derived the axialvector-pseudoscalar 
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and axialvector-axialvector correlators. It turned out that these results not only give a 
smooth interpolation between the e and p regimes but also keep the axial Ward-Takahashi 
identities at an arbitrary choice of quark masses. It will be important to check if this 
simplified prescription is valid for the other quantities like three or four point functions. 
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Appendix A: ^ correlators in finite volume 

Integrals over ^ fields are expressed by A(x, M^) and G{x, Mf, M|) defined by Eqs. (fTTj) 
and f[T^ . Here we summarize useful formulas in the calculation of meson correlators. 

We first note that even simple (three-dimensional) integrals and derivatives of them have 
unusual forms like 

/ d'x A(x M^) - ^Q^MM(t-T/2)) _ J_ 
J ^ ^^^^'^^ ) - 2Msinh(MT/2) M^T' 

1 

V'' 



<9jA(x,M^) = M^A{x,M^) + -, (A2) 



due to absence of the zero-mode. 

:(i) 



For the 0{S\ ) contribution, we need 

1 tpx 

E- 

^^{Aix,M',)-A{x,M',)). (A3) 



r _ _ 1 pipx 



Mi - Mi 

which becomes — 9jv/2 A(x, M^) |^/2=^2 in the limit M^ = M|. 
In the same way. 



J dSA{x - y, Ml)G{y, Ml Ml) = ]_ (G(x, M^, M|) - G(x, M^, M|)) , (A4) 
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which can be expressed in two different ways: 

^—2 {G{x, Ml Ml) - G{x, Ml M, 



M| - M| 

^2 (^(^' Ml Ml) - G{x, Ml, Ml 



3) ; ' 



For contributions, we use 

/ A d'z A{x - y, MlA{y - z, M|) A(^, M^) = 

(^A{x,Ml) - A{x,Ml 



(p2 + M2)2(p2 + M|) 



(M| - Ml)' 



M| - 



dApA{x,M^ 



M'^=Ml 



whose degenerate hmit, Mf = M|, becomes {d^pY A{x, M'^)\j^,[2^j^f2. We also need 

j d'y d^zA^x - y, MlG{y - z, Ml M|) A(z, Ml 

G{x, Ml Ml + M|, M|) - G{x, Ml M|) - Mf, M|) 



(M|-M2)(M|-M2) 



which becomes in the hmit M| = M|, 



_ ^(x, M^, Ml + M^, Ma') - 2G(x, M/, M^) 
" (Af| - M2)2 • 

For the disconnected part, we compute 

d'xA{x,MlA{x,Ml) = ^^[^^ {A{0,Ml-A{0,Mi 



2) j ' 



and 



d'^xAix, MlGix, Ml Ml 



^ (G(0, M|) - G(0, M|) 



1 



g(o,m2,m|)-g(o,m2,m2; 



of which divergent part is treated with the dimensional regularization as usual 



(A5) 



(A6) 



(AT) 



(A8) 



(A9) 



,(A10) 



Appendix B: ^ contraction in the pseudoscalar correlator 

Here we summarize the ^ contractions in (P(x)P(0))°°, (P(x)P(0))^°, (P(x)P(0))2'^ and 
(P(x)P(0))Oi. 
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The first leading contribution is given by 



(P(x)P(0)r = - ' ^ (([^0]l2 - [^0^]2l)([^0]21 - [UlU 



4 



Uo 



{AZf, - AZg) (([f/o]i2 - [U^ohi? - {[Uohi - [Ul]uy)uo 

+^E^(^'^a)([f^o]2at/o^k)c/o 

-^G(x,M2,,M22)(([f/o]i2 + [Ulhi){[Uohi + [f/o^]i2))% 
+^A(a;,M,22)(([f/o]i2r + ([f/J]2inao 

S2 



G(x,M|2,M|2)(([f/o]l2 + [f/o^]2l)')c/o 



4i?2 

where we have used 



4F2 

G(a:,M2,,M2,)(([[/o]2i + [f/o^] 12)^)^/0, (Bl) 



El, (1 - AZg + l^M^) (1 - AZf, + l^M^] = 
Next we calculate the 0{Sj^'') contribution. In this NLO part, we can set Z^-' = Z^j} 



7«J 



1. Note that ^ contractions have to be all connected since the self-contraction is not 
allowed in the NSC vertex in Sj^\ 

Using a notation given in Eq. ( 1721) and the integration formulas given in Appendix El we 
obtain 

S2 . / S 



{Pi^)mr = ^2 ([7^]n + [7^]22)^„ 9M.jA(x,M^ 



2^ 



([7^],l[^o]l.■ + [7^]l,[^oL■l)^o /^(^^ ^2 ) _ ^(^^ ^2 . 

2F^ nij — nil 
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([7^]l2[f/o^]l2 + [f/o]2l[7^]2l) 



2F2 



mi — 7712 



A(x,M22) - A(x,M- 



11, 



([7^]2l[^0121 + [^0]l2[7^]l2)^0 ^2 . _ ^2 

([7^]l2([^0]21 + [^ll]l2) + ([^0]21 + [Ul]l2)[n]2l)u, 

2F'^ mi — m2 

X Mi\, M22) - G{x, M^,, Ml\ 

([7^]2l([^o]l2 + [^o^]2i) + ([^o]i2 + [ul]2i)[nu)u, 

2F'^ 7712 " f^l 



(B2) 

For the O (^{Sj^^y) contribution we have both connected and disconnected parts. Note 

it we can set Z^-' = Zp = Z^li = 1 here, too. 

The connected part (noted by the subscript "con" ) is given by 



(P(x)P(O)) 



20 

con 



2F2 



(['^]lj[^]jl)t/o _^ 



F2 



9a/2 A(x, 



-E '£'"'^'t° (A(x.M;,)-A(x.M: 



.^2 (mi-ma)^ 

(([7^]l2)^ + ([7^]2l)^)^o 
2(m2 — mi)2 



+ 



(A(x, Mi\) + A(x, M2^2) - 2A(x, M- 
^A{x,Ml^,Ml 



(2[7^]l2[7^]2l + ([7^]l2)^ + ([7^]2l)^) 

K-m2)2 ^"^^ 



(B3) 



For the disconnected contribution, we first calculate 



8F2 tl" M2. - M2. V 



AZ 



(B4) 



using Eqs. (lA9p and (lAlOp in Appendix [XI Then we obtain (noted by the subscript "rfis" ] 

1 



(i^(^)i^(0))^°. = ;.[(«(f/o)((5fV)c)^„-(«(f/o))%((5fV)e 



2 
1 

+2 



mUo){{Sf^f)du. - mUo))uASf^f)i,u^ H^,Ml,l (B5) 
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where 



a{Uo) ^ -^[([f/o]i2-[f/ol2i)([t/o]2i-[t/J]i2) 



/3(f/o 



2F2 



11 



-2([^o]i2 

0j22)([f^]22 - 



Ojlly 



(B6) 
(B7) 



Since AZj^ rapidly decreases as the mass reaches the e regime, the contribution is 
important only deeply inside the p regime. Therefore, we can perturbatively perform this 
part of the Uq integral in advance. Using the technique presented in Appendix |Dl the 
calculation is given by 

(a(f/o)7^,,7^,,)^„ - (a(f/o))c/o(7^i,7^,^)^, = 4S'(mi-m2)^ ^^^^^^^ ^ ^^^^^^^ ^ 0(/)j(B8) 



(/3(f/o)7^,,%)c;,,-(/3(f/o))c;„(7^,,7^,,)c/o = Oip'"^), 
where /ij = rriiTiV and we obtain 



(B9) 



where we have used 



22 



(/il+/i2)2 \ 



Uo 



22 



(BIO) 



(Bll) 



Uo 



\ 2 2 
for the later convenience. 

Finally let us calculate the 0{Sj ) contribution. As in the calculation above, using the 
technique in Appendix |Dl we obtain 

/[ ^0 + ^j]ll [^0 + ^0^]22 \ 

2 2 /^^ 



(P(x)P(O)) 



01 



(/il + /i2)2 



X 



12) + f^2 (-AZ. 



2^2 + -J^M!, 



(B12) 



Here we note 



^ / [^0 + [/t],, [^0 + ^J]22 \ 

(/^i + /^2)\ 2 2 /^^ 

In order to obtain the final expression in Eq. fl90|) . we use 

S2 /S 



^cff 



/ 7I2 /7l2\2 



(B13) 



y2 

/ /7I2 /7l2\4/70a I 
4 /il + /i2 V 

neglecting the higher order contributions. 



ccfF I ccff 
A^l + /^2 



(B14) 
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Appendix C: Uq integrals 



The zero- mode Uq integrals of various matrix elements have been calculated in Ref. 16 1. 
Here we summarize the results in our notation for this paper. 



- (^[Uolvv - [Uo]vv) 



Uo 



Uo 



2"' 



2\ ±1 



J {(\.Uo]viV2 ± [f^d]l>2«l) ([^0]l>2t'l i [f4]l'lf2)) 



(CI) 
(C2) 
(C3) 

(C4) 



1 



99 ^^V2^Vi fJ'Vi<Sv2] 
^0 - K2 



(C5) 



Here it is useful to define 



SiSj = lim — — Sj, 



(C6) 



or more explicitly, 



SiSj = < 



lim 



d d 



92 



(C7) 



Note that the partial quenching is performed after the differentiation. Then D's can be 
expressed as 



(C8) 
(C9) 



We note 



m,{S, - 1) ~ 0{p^), 



(CIO) 

(Cll) 



which is useful to simplify our results. 
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We also note that (or I)i2 in the degenerate case mi = m2 = m„) can be written in a 
simpler form than the original definition. Introducing simplified notations for the zero-mode 
partition functions: 



2o = ZQ]^^,({flsea}), 



/^f)2; f^v2, {/^sea}); 



and noting that these partition functions satisfy 



it is easy to show 



for any i. We then obtain 

V = 



lim = Zo, 

lim ^2(m,/ib2|Ai^l,/ii;2) = 

lim Z2{fibl, f^b2\fl'vl, fJ'v2) = 2l(/^fel|/^i;2)5 



+ — I 2,2{fJ'bl, f^b2\f^vl, fJ'v2} 



J__d d_ 



which is used to obtain expressions in Eqs. (I124p and Eq. ( ]13ip . 



(C12) 

(C13) 
(C14) 



(C15) 
(C16) 
(C17) 

(C18) 



(C19) 



Appendix D: Uq integrals in the p regime 

In our calculation, we sometimes encounters a situation that the zero-mode integrals are 
needed only in the perturbative p regime. It is not impossible to nonperturbatively perform 
the zero-mode integrals even in such cases, but it is more convenient to go back to the 
perturbative analysis to obtain the final results in a simple form. 

Let us start with an expansion of the Uq field: 

Uo = exp h-^ ) " + ~ + ■ ■ ■ ' ^ ^ 



and give a Feynman rule for 







{MiA^ohi) = Sii6jkj^. (D2) 



45 



Note that it reproduces the ordinary propagator in the p expansion together with A(x, Mfj). 
It is here important to note that .^o is an element not of SU{N) but of U{N) Lie algebra 
and there is no diagonal contribution like non-zero mode ^ has^. Then we can calculate the 
zero-mode integrals in the p regime as 

{[UoUUo]ki)uo = -SaSjk^^ + 0{p'), (D3) 
{[UoUU^oUuo = +SuSjk^^ + 0{p'), (D4) 



(D5) 

{[UoUulUUo + ulu)^ {[u.uulh)^ {[u, + ulu)^ / 1 Y 

~2 ' 1 ' = ^ ^'"^ [t^J ■ 

(D6) 

These results can be, of course, confirmed by directly performing the exact group integrals 
and then taking the asymptotic expansion in large mjSys. 



Appendix E: Axialvector-Pseudoscalar correlator in the pure e regime 

In this appendix we present the axialvector-pseudoscalar correlator in the e regime, which 
is, to our knowledge, not found in the literature. 

Since ~ C^(c^) is deep inside the e regime, we can neglect the meson mass in the Z 
factors: let us remove the superscripts and use notations such as Zm, Zp. We also note 
Seff = ^Zl^Z\ and AZfa = to NLO in the e regime. 

The source terms are then simplified as 

P'\x) = ([f/o]i2 - PlU) 

E U^) ([f^o]i^5,2 + UUlU) Z^ZpiZM? 

^ This argument is subtle for the summation over topology whose NNLO contribution produces {Q"^) / ^ifij = 
^/ fJ'ifJ'jiJ^f 1//^/): which comes from the diagonal contribution, {[£.o]ii[^o]jj)- Fortunately, however, only 
off-diagonal contributions are needed in the calculation of this paper, and we can therefore ignore this 
subtlety. 
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E[e'(^)]f'' {[Uo]uS,2 - SMh) , (El) 
id 

P'^\x) = (1 o 2), (E2) 
and the axialvector sources can be similarly written as 



i 
2 



121 



= (1^2). (E4) 

Note that the mass term is now an NLO contribution, which can be treated as a pertur- 
bative interaction term and one can omit the mass in the Feynman rule for ^: 

{^ijix)^M{y))^ = 5ii5jkA{x -y,0) - Sij5kiG{x - y, 0, 0), (E5) 

We therefore replace Sj^^ by 

Si^ d'xTr [n'lei^r'''] , (E6) 



where 



n' = M^Uo + U^oM + ^1. (E7) 

Since Sj ~ 0{e^), it is sufficient to calculate 

(Ao(x)P(O)) = i [{Al\x)P'\0) + Al\x)P''{0)r + {Al'{x)P'\0) + Al'{x)P''{0)r 

+ (1 ^ 2). (E8) 

Noting ([9oa - e5oe]?^^(x)[e']^^^(0))^ = 0, and (see Ref. [l 

{[UoMUo]n)uo = ^1 - ^(A^/ + Q){[Uo]n)uo, (E9) 
([f/o^A^f/o^]n)f/o = mi - ^(iV^ - Q){[Ul]n)uo, (ElO) 

ti A f\TT 1 I rrrti 



= i (2 ([t/o]l2 + [U^hi) ([Uohi + [U^oU) + 2 ([f/o]l2 - [t/^]2i) ([Uohl - [t/J]l2 



+ ([t/o]i2 + [U^oWY + + [U^oUy - {[Uo]i2 - [U^ohif - {[Uohl - [U^oW 



(Ell) 
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and using the integration formulas in Appendix |X1 we obtain the correlator, 



{A,{x)pm 



^cff 



'^f + 1 ^0 A(x, Ml,) 



/il /i2 



cff 



Sf + S: 



cff 
2 



1 + !?■ 



cff 

" /if /if . 



5oA(x, 0) 



M^^S^^^ — ^9m29o(G(x, 0) + 0, M^)) 

/il - /i2 



(E12) 



M=0 



Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961). 
S. Weinberg, Phys. Rev. 166, 1568 (1968). 

J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984); Nucl. Phys. B 250, 465 (1985). 
J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987). 
J. Gasser and H. Leutwyler, Nucl. Phys. B 307, 763 (1988). 

C. Bernard [MILC Collaboration], Phys. Rev. D 65, 054031 (2002) |arX iv:hep-lat/0111051| . 



[1] 
[2] 
[3 
K 
[5] 
[6] 

[7] 
[8] 
[9 
[10 

[11 
[12 
[13 

[14] P. H. Damgaard, P. Hernandez, K. Jansen, M. Laine and L. Lellouch, Nucl. Phys. B 656, 226 

(2003) |arXiv:hep-lat/0211020| . 
[15] P. Hernandez and M. Laine, JHEP 0301, 063 (2003) |arXiv:hep-lat /0212014| . 
[16] P. H. Damgaard and H. Fukaya, Nucl. Phys. B 793, 160 (2008) |arXiv:0707.3740l [hepdat]]. 
[17] F. Bernardoni and P. Hernandez, JHEP 0710, 033 (2007) [a rXiv:0707.3887] [hepdat]]. 
[18] F. Bernardoni, P. H. Damgaard, H. Fukaya and P. Hernandez, JHEP 0810, 008 (2008) 

|arXiv:0808.1986] [hepdat]] . 
[19] G. Akemann, F. Basile and L. Lellouch, JHEP 0812, 069 (2008) [arXiv:0804.3809] [hepdat]]. 

48 



G. Colangelo, S. Durr and C. Haefeh, Nucl. Phys. B 721, 136 (2005) ^arXivihepdat /0503014] . 
J. Gasser and H. Leutwyler, Phys. Lett. B 188, 477 (1987). 

H. Neuberger, Phys. Rev. Lett. 60 (1988) 889. 

F. C. Hansen, Nucl. Phys. B 345, 685 (1990); F. C. Hansen and H. Leutwyler, Nucl. Phys. B 
350, 201 (1991). 

P. Hasenfratz and H. Leutwyler, Nucl. Phys. B 343, 241 (1990). 
H. Leutwyler and A. Smilga, Phys. Rev. D 46, 5607 (1992). 

P. H. Damgaard, M. C. Diamantini, P. Hernandez and K. Jansen, Nucl. Phys. B 629 (2002) 
445 far Xiv:hepdat/01 12016] . 



[20] A. Shindler, Phys. Lett. B 672, 82 (2009) |arXiv:0812 .22511 [hep-lat]]. 
[21] O. Bar, S. Necco and S. Schaefer, JHEP 0903, 006 (2009) [arXiv: 08 1 2 . 2403) [hep-lat] ] . 
[22] F. Bernardoni, P. Hernandez and S. Necco, JHEP 1001, 070 (2010) |arXiv:0910.2537l [hep- 
lat]]. 

[23] O. Bar, S. Necco and A. Shindler, JHEP 1004, 053 (2010) [ arXiv:1002. 15821 [hep-lat]]. 
[24] F. Bernardoni, N. Garron, P. Hernandez, S. Necco and C. Pena. [arXiv:1008.1870l [hep-lat]. 
[25] C. Lehner, S. Hashimoto and T. Wettig, JHEP 1006, 028 (2010) [arXiv: 1004.5584) [hep-lat]]. 
[26] T. DeGrand, Z. Liu and S. Schaefer, Phys. Rev. D 74, 094504 (2006) [Erratum-ibid. D 74, 

099904 (2006)] [ arXiv:hep-lat/0608019] . 
[27] C. B. Lang, P. Majumdar and W. Ortner, Phys. Lett. B 649, 225 (2007) 

|arXiv:hep-lat/06110l0l . 

[28] P. Hasenfratz, D. Hierl, V. Maillart, F. Niedermayer, A. Schafer, C. Weiermann and M. Wein- 

gart, JHEP 0911, 100 (2009) |arXiv:0707.007n [hep-lat]]. 
[29] H. Fukaya et al. [JLQCD collaboration], Phys. Rev. D 77, 074503 (2008) [ arXiv:0711.4965] 

[hep-lat]]. 

[30] A. Hasenfratz, R. Hoffmann and S. Schaefer, Phys. Rev. D 78, 014515 (2008) [arXiv: 0805. 2369) 
[hep-lat]]. 

[31] A. Hasenfratz, R. Hoffmann and S. Schaefer, Phys. Rev. D 78, 054511 (2008) [arXiv:0806.4586l 
[hep-lat]]. 

[32] O. Bar, S. Necco and S. Schaefer, PoS LAT2009, 078 (2009) |arXiv:0910.2372l [hep-lat]]. 

[33] K. Jansen and A. Shindler, PoS LAT2009, 070 (2009) |arXiv:0911.193n [hep-lat]]. 

[34] P. H. Damgaard and H. Fukaya, JHEP 0901, 052 (2009) |arXiv:0812.2797] [hep-lat]]. 

[35] P. H. Damgaard and S. M. Nishigaki, Nucl. Phys. B 518, 495 (1998) |arXiv:hep-th/9711023] . 

[36] T. Wilke, T. Guhr and T. Wettig, Phys. Rev. D 57, 6486 (1998) |arXiv:hep-th/9711057| . 

[37] G. Akemann and P. H. Damgaard, Nucl. Phys. B 528, 411 (1998) |arXiv:hep-th/9801133| . 

[38] A. V. Smilga and J. Stern, Phys. Lett. B 318, 531 (1993). 

[39] J. C. Osborn, D. Toublan and J. J. M. Verbaarschot, Nucl. Phys. B 540, 317 (1999) 
|arXiv:hep-th/9806110] . 

[40] H. Fukaya, S. Aoki, S. Hashimoto, T. Kaneko, J. Noaki, T. Onogi and N. Yamada [JLQCD 
collaboration], Phys. Rev. Lett. 104, 122002 (2010) [Erratum-ibid. 105, 159901 (2010)] 
|arXiv:0911.5"555l [hep-lat]]. 

49 



[41] H. Fukaya et al. [JLQCD and TWQCD collaborations], [arXmT0T23052] [hep-lat]. 

[42] P. H. Damgaard and K. Splittorff, Nucl. Phys. B 572, 478 (2000) |arXiv:hep-th/99121 46]. 

[43] P. H. Damgaard and K. Splittorff, Phys. Rev. D 62, 054509 (2000) |arXiv:hep-lat70003017] . 

[44] P. H. Damgaard, Phys. Lett. B 476, 465 (2000) [arXiv:hep-lat /0001002] . 

[45] J. Wess and B. Zumino, Phys. Lett. B 37, 95 (1971). 

[46] E. Witten, Nucl. Phys. B 223, 422 (1983). 

[47] K. Splittorff and J. J. M. Verbaarschot, Phys. Rev. Lett. 90, 041601 (2003) 

|arXiv:cond-mat/020959 4 1 . 
[48] Y. V. Fyodorov and G. Akemann, JETP Lett. 77, 438 (2003) [Pisma Zh. Eksp. Teor. Fiz. 77, 

513 (2003)] |arXiv:cond-mat/021 0647]. 
[49] K. Splittorff and J. J. M. Verbaarschot, Nucl. Phys. B 683, 467 (2004) |arXiv:hep-th/0310271] . 
[50] J. Lenaghan and T. Wilke, Nucl. Phys. B 624, 253 (2002) [arXiv:hep-th/0108166] . 
[51] S. Aoki and H. Fukaya, Phys. Rev. D 81, 034022 (2010) [arXiv: 0906. 4852] [heo-lat]]. 
[52] Y. Y. Mao and T. W. Chiu [TWQCD Collaboration], Phys. Rev. D 80, 034502 (2009) 

|arXiv:0903.2"T46] [hep-lat]]. 
[53] S. R. Sharpe and N. Shoresh, Phys. Rev. D 62, 094503 (2000) [arXiv:hep-lat/0006017] . 
[54] J. Bijnens, N. Danielsson and T. A. Lahde, Phys. Rev. D 73, 074509 (2006) 

|arXiv:hep-lat/0602003] . 
[55] M. Luscher, Commun. Math. Phys. 104, 177 (1986). 
[56] H. Leutwyler, Phys. Lett. B 189, 197 (1987). 

[57] P. Hasenfratz and F. Niedermayer, Z. Phys. B 92 (1993) 91 |arXiv:hep-lat/9212022| . 
[58] P. Hasenfratz, Nucl. Phys. B 828, 201 (2010) |arXiv: 0909. 34191 [hep-th]]. 
[59] W. Bietenholz et al, Phys. Lett. B 687, 410 (2010) |arXiv:1002.1696l [hep-lat]]. 
[60] M. Weingart. larXiv:1006.5076l [hep-lat]. 

[61] F. Niedermayer and C. Weiermann, Nucl. Phys. B 842, 248 (2011) [arXiv: 1006.5855) [hep-lat]] . 
[62] W. Bietenholz et al. [QCDSF Collaboration], J. Phys. Conf. Ser. 287, 012016 (2011) 
|arXiv:1103.3"3TT] [hep-lat]] . 



50 



